Surface Phase Diagrams for Wetting on Heterogenous Substrates 



C. Rascon* and A. 0. Parry 
Department of Mathematics, Imperial College 
180 Queen's Gate, London SW7 2BZ, United Kingdom 



o 
o 



> 

oo 



o 

O 



o 
o 



We propose a simplified description of fluid adsorption 
on heterogenenous micropatterned substrates. Using this ap- 
proach, we are able to rederive results obtained earlier us- 
ing effective interfacial Hamiltonian methods and predict a 
number of new examples of surface phase behaviour for both 
singly and periodically striped substrates. In particular, we 
show that, for a singly striped system, the manner in which 
the locus of surface unbending phase transitions approaches 
the pre-wetting line of the infinite pure system, in the limit 
of large stripe widths, is non-trivial and sensitive to several 
characteristic lengthscales and competing free-energies. For 
periodic substrates, we investigate finite-size deviations from 
Cassie's law for the wetting temperature of the heterogeneous 
system when the domain sizes are mesoscopic. 



PACS numbers: 



;.08.B, 68.43, 68.35.R, 68.35.M 



I. INTRODUCTION 



Fluid interfacial phenomena on flat, homogeneous sub- 
strates have drawn enormous attention over the last few 
decades [Q. However, present experimental methods al- 
low an ever increasing control over the shape and chemi- 
cal composition of solid surfaces and pose more demand- 
ing challenges for both theory and experiments. These 
advances make the study of adsorption on non-planar and 
heterogeneous substrates inevitable. Recent experimen- 
tal and theoretical studies have shown that structured 
substrates exhibit a variety of novel adsorption proper- 
ties which not only promise to be of importance to future 
technologies such as microfluidics Q , but are also of fun- 
damental interest to statistical physics ||. In particular, 
the effect of chemical inhomogeneities has been addressed 
recently in different contexts, for example: contact angles 
of liquid drops , droplet spreading jjj , morphological 
phase transitions ||, three phase contact line ||, Cassie's 
law drop shapes |l2|| , construction of magnetic 

materials [ fL3| , microscopic packing Q), liquid channels 
fL5| , polymer blends p6| , ^7t , and dewetting |Q, among 
others. 

In this paper, we concentrate on the equilibrium wet- 
ting properties of flat but chemically structured sub- 
strates. Specifically, we study liquid adsorption in two 
representative micropatterned systems: a homogeneous 
substrate with a single distinct chemical stripe and a sub- 
strate comprising parallel stripes of two different materi- 
als. These systems have been partially studied by other 
authors with the help of density functional theory [[t4 19 
and effective interfacial methods I2CH23]. Nonetheless, 



the intrinsic complexity of those formulations does not 
allow a thorough exploration of the vast space of possible 
surface phase behaviour these systems offer. We propose 
here a minimal model to understand the general wetting 
properties of these systems. The model, described in 
section II, appears to capture the essential physics of the 
problem in a qualitative manner and recovers former re- 
sults obtained with more sophisticated methods for par- 
ticular cases. The great advantage of the minimal model 
is that it allows us to obtain (sections of) global surface 
phase diagrams and explore what kind of surface phase 
behaviour is possible in such heterogeneous systems. As 
we shall show, in spite of the simplicity of the model, 
the ensuing surface phase behaviour is extremely rich. 
In section III, the model is applied to a substrate with a 
single stripe of a different material, while in section IV, 
a periodic substrate made of parallel stripes of two ma- 
terials is studied. We finish our article with a summary 
of our main conclusion and a discussion of further work. 



II. MINIMAL MODEL 

Consider two flat substrates made of different mate- 
rials (labelled 1 and 2) whose adsorption properties can 
be described, when homogeneous, by effective potentials 
W\(l) and respectively [p3[ . When each of these 

substrates is exposed to a near saturated g elS, el liquid 
layer is adsorbed on its surface, whose thickness t\ (or 
£2) is given by the position of the absolute minimum 
of the corresponding effective potential (the superscript 
7r denotes behaviour characteristic of the pure homoge- 
neous substrate). We suppose the wetting temperatures 

of each of these systems are T„ and \ respectively. 
Consider now a flat substrate A consisting of domains 
Ai and A2 (not necessarily connected) made of materi- 
als 1 and 2, respectively. When this microprinted sub- 
strate is exposed to a gas, a liquid layer is adsorbed on its 
surface, whose equilibrium thickness profile l{x, y) must 
show considerable spatial dependence (especially strong 
close to the boundaries between the domains) . For meso- 
scopic domains, this dependence can be captured by an 
effective interfacial description: the equilibrium profile 
of the adsorbed liquid layer will minimise the following 
free-energy functional 
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where £ is the liquid-vapour interfacial tension. This 
model assumes a sharp crossover in the local form of the 
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effective potential at the domains boundaries and is it- 
self an approximation to a more realistic description in 
which this crossover is smooth ||(J. However, this as- 
sumption certainly does not alter the global structure of 
the surface phase behaviour ||. To continue, and for the 
sake of simplicity, we only consider systems translation- 
ally invariant along the y axis (see Fig. 1). Therefore, 
the equilibrium profile will follow from the minimisation 
of the somewhat simpler functional 

H[l] = [dx | l 2 + [ dx Wt[l) + [ dx W 2 (£), (2) 
J c 2 J Cl J C2 

where the dot represents differentiation with respect to 
x, and C, C\ and C 2 are the appropriate linear domains 
corresponding to sections of A, Ai and A2, respectively. 
This approach has been employed by a number of au- 
thors to obtain numerical solutions of the equilibrium 
profile for a substrate with a single stripe of a different 
material (Fig. 1(a)) @-^2|. Their results show that if 
the width of the stripe L\ is very large, the thickness 
of the adsorbed layer is essentially the same as in the 
pure systems within each of the domains {P[ in Ai , £ 2 m 
A2), and varies abruptly in the domain borders. In other 
words, for large stripe widths, the minimisation of (Q) is 
achieved by a local minimisation of the second and third 
terms of the functional. The (positive) contribution of 
the first term in (||), which originates in the domain bor- 
ders, is practically negligible. We can estimate this latter 
contribution by noticing that the characteristic length for 
lateral variations of the interface profile is given roughly 
by the transverse correlation lengths £,| = y/T,/W"(£j) 
of the homogeneous substrates i = l,2. Thus, for a single 
step, we can write 

Jdxp 2 « \ (3) 

with k = )■ For smaller values of L\, how- 

ever, this energetic contribution becomes comparable to 
that arising from the second term of (||) and the inter- 
face adopts distinct configurations in order to minimise 
the free-energy of the system. This gives rise to certain 
surface phase transitions, denoted surface condensation 
or unbending transitions, in which different liquid interfa- 
cial configurations coexist [pp|- [2^j2^1 . In order to calcu- 
late the loci of these transitions in the T-fx diagram, the 
free-energy functional (||) must be minimised for different 
values of the temperature T and the chemical potential ji. 
Although this has been done for a specific case pjj ], the 
vast space of parameters involved and the arduous task of 
minimising the functional (g) does not allow a thorough 
exploration of the physics of the system and impedes a 
systematic study of more complex systems, for instance, 
a periodic array of stripes (see Fig. 1(b)). 

Here we propose a minimal model to describe the over- 
all phenomenology of the above mentioned surface un- 
bending transitions in a simple and illuminating manner. 



To do that, we approximate the free-energy functional 
by a free-energy function of a finite number of vari- 
ables. By doing that, we lose the detailed description 
of the shape of the adsorbed layer (which follows from 
the minimisation of Eq. @), but, as explained here, we 
capture the fundamental physical mechanism for the sur- 
face phase transitions underlying (Q). To start, we define 
a collective coordinate li for every connected region of 
each of the domains. This coordinates will account for 
the average thickness (in a loose manner) of the adsorbed 
layer within that connected region and it is defined im- 
plicitly as follows 

LiW a {£i) = J dx W a (l), (4) 

Li being the width of the connected region and a = 1 
or 2, depending on the domain that region belongs to. 
Within the minimal model, we consider these collective 
coordinates as independent entities, to be varied in or- 
der to minimise the free-energy. This accounts for the 
second and third integrals of (El) which are replaced by a 
sum over terms Li W a (£i), one for each connected region. 
To obtain a full description in terms of these variables, 
we simply estimate the contribution of the first term of 
(||) similar to that shown in Eq. (||) but using the col- 
lective coordinates. Thus, the free-energy functional is 
reduced to a function of a small number of real variables. 
Specifically, for the periodic array of stripes (Fig. 1(b)), 
we write the free-energy per period as 

H(£ u e 2 ) = k (4-4) 2 + L x W x {h) + L 2 W 2 (t 2 ). (5) 

In accordance with our aim of developing a minimal de- 
scription, we consider k constant throughout the paper 
(independent of the temperature and the chemical po- 
tential) . 

The free-energy for a substrate with a single stripe 
(Fig. 1(a)) can be obtained from (|^) in the limit L 2 ~^oo. 
In this case, the last term of (^) becomes the most im- 
portant and it must be minimised fully, yielding £ 2 =£ 2 . 
Substracting this (infinite) constant contribution from 
the free-energy, the average thickness of the adsorbed 
layer on the stripe will follow from the minimisation of 

H(h) =K(h-t 2 T ) 2 +LiW 1 {h). (6) 

The generalization of this procedure to more sofisti- 
cated geometries is straightforward. As an example, the 
average thickness l\ of the adsorbed layer on a circle of 
radius R made of material 1 on a substrate made of ma- 
terial 2 could be obtained from the minimization of 

H{t x ) = nRn^-q) 2 + ttR 2 Wy{h), (7) 

where an appropriate k now accounts for the energy due 
to the curvature of the interface. It is likely that this 
procedure will fail to work for certain intricate shape 
domains where more complicated behaviour may be ex- 
pected. 



2 



We stress here that this approach does not intend to 
be a quantitative description of these phenomena but to 
capture the essential features of heterogeneous adsorp- 
tion, which, due to the complexity of the available mod- 
els, remain mainly unexplored. 

III. SINGLE STRIPE 

A. General Considerations 

Here we apply the model described above to a homo- 
geneous flat substrate with an infinite length stripe of 
width L\ (Fig. 1(a)). Before presenting the results, how- 
ever, we make a number of pertinent remarks concerning 
general aspects of the surface phase behaviour: 

• For small values of the stripe width, corresponding 
to the limit L\ — > 0, we expect that the thickness 
of the layer adsorbed on the stripe will tend to the 
thickness of the liquid layer absorbed on the outer 
substrate, i.e. 

• If the stripe width is large, corresponding to the 
limit L\ — > oo, the thickness of the layer adsorbed 
on the stripe will tend to the thickness of the liquid 
layer absorbed on a homogeneous flat substrate of 
material 1, i.e. 

• If the thickness of the adsorbed layers on both in- 
finite substrates is the same (l\ = £J), the func- 
tional (||) is fully minimised by a flat solution 
l{x) — l\ — t\. The loci of the points for which 
l\ = constitute a line in the T — p phase di- 
agram and separates interfacial profiles with dif- 
ferent convexity |2q| . Observe that the flat solu- 
tion corresponds to the minimisation of each term 
of the functional and cannot be improved by any 
other configuration. It follows that such a solution 
will be the global equilibrium for arbitrary values 
of the stripe width L\. In turn, this implies that 
no first-order surface coexistence line can cross the 
line defined implicitly by P[ = i^- As explained 
below, this simple observation severely constrains 
the possible topology of the surface phase diagram. 
Intriguingly, the loci of points satisfying l\ — l\ 
represents surface phases which are fully homoge- 
neous on the heterogeneous substrate. 

• Any first-order phase boundary corresponding to 
a line of coexistence in the T — p plane satisfies a 
Clausius-Clapeyron-like equation ||26[ : 

(dp\ _ AS 

\ / COGX 

where AS and Ar are the differences in surface en- 
tropy and adsoption between the coexisting phases, 



respectively, which follow in the usual manner from 
the free-energy of the system u)(T, p) = mini? as 
S=-(dLu/dT)^ and V=-{du)/dp) T . 

For our purposes, it is more convenient to measure 
the slope of this line with respect to that of the 
liquid- vapour coexistence curve p LV (T). Thus, we 
rewrite Eq. (||) as 

= (9) 

where Ap = (p LV — /z), AT = T — T w (for a given 
wetting temperature), S = — (dcu/d AT)^ and T = 
-(<9w/5A/x) A T = -r. 

B. Minimal Model 

The equilibium value of the average thickness 
£i(T, p, L\) follows from the minimisation of the mini- 
mal free-energy (@), yielding 

2K(l x -tl) + L x W[(l x ) =0. (10) 

This equation has a straightforward graphical interpreta- 
tion from which one can elegantly glean the possible sur- 
face phase behaviour. To see this, we rewrite the equation 

as 

W[(h) = ^ {ll-h), (11) 

corresponding to the intersection of a straight line with 
a curve. The general form of the graphical construction 
is illustrated in Fig. 2 for a specific choice of effective 
potential. Note the following features: 

• The model always shows at least one solution. Mul- 
tiple solutions correspond to the presence of first- 
order surface phase transitions. 

• The straight line has a negative slope inversely pro- 
portional to the stripe width and necessarily crosses 
the horizontal axis at 

• The curve always intersects the horizontal axis at 
i\ (but may exhibit other crossings) and tends to 
A/i Sp as l\ — > oo. Here, bp represents the difference 
in bulk densities of the coexisting liquid and gas 
states. 

• For any given temperature and chemical potential, 
only three scenarios are possible within the graph- 
ical construction. These correspond to the straight 
line intersecting the curve at one, two or three 
points. In the first scenario [Fig. 2(a) or (b)], the 
intersection constitutes the unique solution of Eq. 
(pr|). In the second case [Fig. 2(c)], the straight line 
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must be tangential to the curve at a point (repre- 
senting a spinodal) and crosses the curve at a sec- 
ond point, which constitutes the stable solution. In 
the third case, the central intersection represents an 
unstable solution whilst the other two are mechan- 
ically stable. Among them, the stable (metastable) 
solution will have the lower (higher) free-energy. If 
both solutions have the same free-energy, two dif- 
ferent thicknesses of adsorbed layer coexist at the 
wall-fluid interface corresponding to a first order 
phase boundary [Fig. 2(d)] (§7|. 

• In the limit L\ — > 0, the straight line is vertical 
and the solution tends to the thickness of the liq- 
uid layer adsorbed on the surrounding substrate, 
ti—ti^, as expected (see previous subsection). Al- 
though the model is not intended to be accurate 
for very narrow stripes, note that it recovers the 
correct result in this limit. 

• In the limit L\ — *■ oo, the straight line is horizon- 
tal and the solution recovers the homogeneous sub- 
strate, t\^ri\ (see previous subsection). 

• If £±=£2, both lines cross at the same point of the 
horizontal axis and there is only one solution inde- 
pendent of the value L\ (see previous subsection) . 

• There is a critical point for a certain combination 
of the parameters of the problem. This point cor- 
responds to the straight line intersecting the curve 
tangcntially at the point of zero curvature (and neg- 
ative slope). See Fig. 2(b). Thus, within this ap- 
proach, a necessary condition for the existence of 
a critical point is the presence of a point in the 
effective potential for which the third derivative is 
zero (the second derivative at that point being neg- 
ative). This condition can be used to determine a 
bound for the critical point (see below). 

Due to the simplicity of the model, the conditions 
for criticality can be written explicitly as 



and 



W["(h) = 0, 



(12) 



(13) 



which, together with Eq. ( |11[ ) determine the critical 
point completely. Observe that Eqs. ([tl]) and Jl2] ) 
constitute together the conditions for the spinodals. 



C. Second-order Wetting Stripe 

To begin we consider the case where the stripe corre- 
sponds to a material that exhibits a continuous (critical) 
wetting transition in the limit of infinite stripe size. No 



assumptions about the outer domain have to be made at 
this moment. Here we concentrate on the most exper- 
imentally relevant case of long-ranged forces and write 
our model effective potential in the usual manner [|l| 



AT 



A 



A/i bp 



(14) 



where we have rescaled the temperature such AT = 
T-T^\ A is a (positive) Hamaker constant and p de- 
termines the range of the forces (p = 2 for Van der Waals 
dispersion forces). 

The effective potential (|lj) models the continuous di- 
vergence of the film thickness as the temperature is in- 
creased towards the wetting temperature at Ap = 0. 
Observe that the third derivative of this potential con- 
sists of only two terms and does not depend explicitly on 
the chemical potential. Recalling the condition for crit- 
icality, Eq. ([II]), it follows that a critical point can only 
occur provided AT < 0. In other words, if this system 
has a critical point, it must occur necessarily at a tem- 
perature below the wetting temperature T, 



(i) 



The conditions for criticality follow immediately from 
Eqs. © to (H): 



ALx{p + l) 



2 k 



1 



ATr 



p + 3 A 

P li n 



AjjL C bp 



2k 



<2 



P 



(15) 



(16) 



(17) 



Figure 3 shows the phase diagram obtained for the ef- 
fective potential (14) and different values of the stripe 
width L\. For the sake of simplicity, the thickness of the 
adsorbed layer in the surrounding substrate, 1%, has been 
considered constant (independent of the temperature and 

the chemical potential). This approximation represents 

(2) 

a system whose wetting temperature T4, is higher than 
the wetting temperature T^P and for which the varia- 
tion in the thickness of the adsorbed layer is negligible 

at temperatures close to 2w . The advantange of doing 
is that one can more easily identify the origin of certain 
features of the surface phase behaviour (see below). The 
microscopic units of length of the problem are given by 
A= (A/k) 1 / 4 . In the calculation, we have fixed the range 
of the forces to be p = 2 and the thickness l\ — 4A. The 
qualitative features of the phase diagram described here 
are not in any way specific to these values. 

For L\ — 12A, the system undergoes an unbending 
transition represented by a continuous thick line ending 
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in a critical point (see Fig. 3). Along the line, the ad- 
sorbed layer on the stripe can have either of two different 
thickness, both of them thinner than the surrounding 
adsorbed layer l\. Note that the unbending line joins 
the liquid-vapour coexistence line (A/i = 0) at a certain 
non-zero angle. This is consistent with Eq. (||) and the 
fact that the difference in coverage Ar ~ A£i between 
the coexisting phases at that point is finite. This transi- 
tion has been predicted recently p2| , although the T—fi 
phase diagram had not been described. Observe that the 
critical temperature of this transition lies below as 
predicted above. As the stripe width is increased, the 
length of the unbending line decreases and for a given 
finite value of the width L\ = L\ w 53A, it vanishes (white 
circle in Fig. 3). This phenomenology is identical to that 
described for a homogeneous but corrugated substrate 
pij , showing that the physics involved is essentially the 
same: the transition takes place due to a balance between 
the free-energy associated with the direct intemolecular 
interaction with the substrate and the energetic cost of 
increasing the area of the liquid-vapour interface. For 
this particular model, Eq. ([l4j), the transition vanishes 
for 



t 



2 k 



A(p+1) 



(p+1)^ 



p+3 



at a temperature 



AT 1 * 



(p + 3) 2 A 
p(p+l) 



»7T 

I. 2 



(18) 



(19) 



The loci of the critical points are represented in the 
same figure by a grey thick line. For L\ > L\, this line 
extends further into the metastable region of the phase 
diagram (dashed grey line) and, in the limit L x — > oo, 
merges asymptotically with the wetting critical point at 
2w (large black circle in Fig. 3). The presence of a 
(metastable) critical point close to the coexistence line 
will affect the adsorption properties in the stable part 
of the phase diagram. For instance, the adsorption 
on the stripe will show a strong increase close to that 
(metastable) critical point. 

The loci of the points for which P[ =l\ are represented, 
in the same figure, by a straight black dotted line. Whilst 
this line does not represent any surface phase transition 
or singular behaviour, it does significantly influence the 
possible form of the phase diagram. To see this, recall, 
as mentioned at the beginning of this section, that the 
adsorbed layer is flat along that line for any value of the 
stripe width L\. Thus, neither the loci of the critical 
points nor any unbending phase boundary can cross the 
line £? = £%. 



D. First-order Wetting Stripe 

We now suppose that the stripe corresponds to a ma- 
terial that exhibits a first-order wetting transition in the 



limit of infinite stripe size. For this case, we use an effec- 
tive potential 



A 2 
4B 



-AT 



where AT = T — \ A and B are (positive) Hamaker 
constants and, as earlier, p determines the range of 
the forces. This potential models the discontinuous di- 
vergence of the film thickness l\ at and an as- 
sociated pre-wetting line off coexistence for AT > 0. 
This pre-wetting line ends in a pre-wetting critical point 
(AT pw , A/i pw ) where 



1 



A 



B 



£p+l £p+2 



+ ApSp 



(20) 



AT = 1 ^ 
pw 2p(p + 3) B 



Ap pw dp = 



2B 



P 



1 



1 A 



1 P+3 



2 (p + 3) B 



At this point, the wetting layer thickness is 

/ = 2 (p + 3) B 
A 



P 



1 



(21) 



(22) 



(23) 



which will be useful in the following analysis. 

A necessary condition for criticality in the finite stripe 
system is that the effective potential has a point where 
the third derivative vanishes, (Eq. (|i~3|)). This provides 
an upper bound in temperature for any possible critical 
point in the system, AT C < 3 A 2 / {<ip {p+4)B). In contrast 
with the previous case (section [II C), a critical point can 
appear in this system at temperatures above the wetting 
temperature T4 (but below a certain threshold). The 
conditions for the critical point, Eqs. (|TT|)-([T^) , read 

2k 



(p + 1) A £ lc - 2 (p + 3) B + — £ 



p+4 



0, (24) 



ATr 



p + 3 A (p + 3)(p + 4) B A 2 



P il r 



p(p + l) £ 



1 



4Ti ' 



and 



2k 

Afi c op= — 



q- p -±le lr 
2 p + i c 



B 



P 



+ 1 e, p+3 ' 



(25) 



(26) 



In the limit L\ — > oo, these equations recover (^l|)- 
(p3|), showing that the critical point of the system must 
merge into the pre-wetting critic al po int for large stripe 
widths, as expected (see section [II A ). Notice that Eq. 
( p4| ) implies that the thickness at the unbending critical 
point l\ c must be lower than the pre-wetting critical 
thickness £ pw for any finite value of L\. This allows us 
to sharpen our upper bound for the finite stripe width 
critical temperature to Tq <T pw for all widths L\. 
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Figures 5-8 show the phase diagrams obtained for the 
effective potentiai ( p0| ) and different values of the stripe 
width L\. Once again, we have chosen A = (A/n) 1 / 4 as 
unit of length and the range of the forces has been fixed 
to p = 2. Besides, due to the large number of parame- 
ters, we have fixed the Hamaker constant B = XA. As 
in the previous calculation, the thickness of the adsorbed 
layer in the surrounding substrate, has been consid- 
ered constant (independent of the temperature and the 
chemical potential). However, in contrast to the case of 
critical wetting, different values of l\ l ea d to quite dif- 
ferent surface phase behaviours. As we shall show, the 
discriminating parameter here is the ratio of £\ to £ pw . 

The phase diagram of this system for £\ = 6A and 
L\ = 2000A is plotted in Fig. 5. As before, we find a 
line of first-order undending phase transitions ending at 
a critical point (black continuous line). Along this line, 
phases with distinct adsorptions coexist (see insets) . The 
line i\ =£2 represents the loci of points for which the ad- 
sorbed layer thickness is homogeneous even though the 
substrate is heterogenous. The line serves to separate 
interfacial configurations with different convexities 
An extension of this line goes beyond the pre-wetting 
line and joins the coexistence line. Along this extension 
(long dashed line), the liquid interface is also flat due to 
the presence of a metastable minimum or even a maxi- 
mum of the effective potential (^0|) for £ — £2 [Hi • Note 
that this scenario differs from that occuring for a second- 
order wetting effective potential. In that case, the line 
£\ = £2 never hits the coexistence line for any value of 
the stripe width (see Fig. 3) and, therefore, coexistence 
always takes place between adsorbed layers thinner than 
the layer adsorbed on the surrounding substrate. In ei- 
ther case, we designate the transition as unbending since 
it originates from the same balance between the energy 
associated to the interaction with the substrate and the 
energetic cost of increasing the area of the liquid- vapour 

interface HII. 

Observe also in Fig. 5 that the unbending line joins the 

liquid- vapour coexistence line (A/i = 0) at a certain non- 
zero angle whilst the pre-wetting line approaches that line 
tangentially. This is due to the fact that the change in 
the order parameter, £\, across the unbending transition 
is finite for the unbending transition while is infinite for 
the pre-wetting transition (because the system is wet for 

T>Tw at liquid- vapour coexistence) [see Eq. ([)])]. 

Figure 6 shows the phase diagram of this system for 
£2 = 10A and different values of the stripe width L\. As 
L\ increases, the coexistence line approaches the pre- 
wetting line. A detailed description of the merging of 
both lines is given below. Notice the presence of the line 
£\ — £2 above the pre-wetting transition and the fact that 
the unbending transition does not cross that line for any 
value of L\. 

The phase diagram for l\ — 4A and different values 
of the stripe width L\ is presented in Fig. 7. As ex- 
pected, the unbending line merges the pre-wetting line 



in the limit L\ — >oo. Note that the merging is not mono- 
tonic: as L\ increases, the unbending line crosses the pre- 
wetting line and merges from above, close to the coexis- 
tence line, and from below, close to the pre-wetting criti- 
cal point (see inset). The locus of the unbending critical 
points is similarly non-monotonic for two reasons. On 
one hand, for small values of Li, the unbending transi- 
tion occurs at (AT, A/i where the potential (^0|) has only 
one minimum. Consequently, the loci of the unbending 
critical points resemble those shown in Fig. 3 for criti- 
cal wetting. On the other hand, for larger values of L\, 
the unbending transition takes place at (AT, A/i) where 
the potential ( ^p| ) shows the characteristic double mini- 
mum of pre-wetting phase coexistence and the loci of un- 
bending critical point re-routes towards the pre-wetting 
critical point. The upshot of this is that pre-wetting only 
has influence on the surface phase diagram for sufficiently 
large stripe widths L\. 

The contrast between these mechanisms is clearly mag- 
nified for thinner adjacent wetting layers. Fig. 8 shows 
the phase diagram obtained for £\ = £ pw = 10A/3. Ob- 
serve that the loci of the unbending critical points now 
cross the liquid-vapour coexistence curve and, for a cer- 
tain range of width stripes (94A < L\ < 938A), the un- 
bending phase transition is not present since it is pre- 
empted by bulk condensation. The dissappearence of the 
unbending transition as the stripe width increases (for 
L\ « 93. 6A) resembles the behaviour for critical wetting 
(see Fig. 3). In contrast, when L\ > 937. 5A, the un- 
bending transition re-enters the phase diagram (at tem- 
peratures above T^) and, in the limit L± — > 00, fuses 
with the pre-wetting line. Notice once more that the be- 
haviour of the unbending line with the stripe width is 
non-monotonic: after re-entering the phase diagram, its 
overall position first rises in temperature and, for larger 
values of L±, lowers towards the pre-wetting line (see 
inset). For the particular effective potential ((2^), and 
within the context of the minimal model, this re-entrant 
behaviour for unbending occurs when 



I 



q < (p + 2) 2 



pw 



(p+l)(p + 3)' 



(27) 



Observe that the disappearence of the unbending tran- 
sition (and further re-entry at higher temperatures) for 
small values of £\ is inevitable if the unbending transi- 
tion is not to cross the iso-adsorption line £\ =£2, since 
this line surrounds the pre-wetting critical point as £2 
decreases (see Fig. 4). For this reason, this feature must 
survive in analysis based on the full interfacial model, 
Eq. (|). 



E. Discussion 

Before we turn to the results for a periodic array of 
stripes, we complete this section with some critical re- 
marks concerning the features of the phase diagrams pre- 
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dieted above. To facilitate this, it will be instructive 
to compare our results with those obtained within the 
context of the full interfacial model (||), by Bauer and 
Dietrich (BD) for a particular stripe system made of a 
material which undergoes a first-order phase transition 
(when infinite) pi} ). Specifically, let us concentrate on 
the T—fi section of the phase diagram shown in Fig. 2 of 
their article. Whilst their diagram resembles our Fig. 7 in 
a qualitative way, there are three features which deserve 
detailed comment: 

• In the limit L\ — > oo, the adsorption characteris- 
tics of the stripe must tend to those of the pure 
substrate, i.e. the unbending line must merge with 
the pre-wetting line. Both approaches capture this 
requirement, although they differ in the manner in 
which the mentioned lines coalesce. The main dif- 
ference lies in the behaviour of the unbending line 
close to liquid- vapour coexistence. In the BD phase 
diagram, the unbending line fuses the pre-wetting 
line monotonically from lower temperatures. In 
contrast, as described above, the minimal model 
predicts a non-monotonic merging: close to liquid- 
vapour coexistence, the unbending line approaches 
the pre-wetting line from higher temperatures, as 
Li^oo (see, for example, the inset of Fig. 7). 

A detailed calculation with the full interfacial 
model (||), included in appendix [a|, shows that ei- 
ther of both situations is possible, and that the 
actual prevalence of one or the other depends on a 
number of factors. However, for the experimentally 
relevant case of dispersion forces, p — 2, we believe 
the scenario described by the minimal model is the 
correct one: the unbending line merges with the 
pre-wetting line from higher temperature close to 
the liquid-vapour coexistence curve. We suspect 
that numerical problems involved in the large-scale 
computation of BD are to blame here. For future 
numerical studies, we note that the behaviour of 
the unbending line close to coexistence is sensitive 
to the range of the forces. Consequently, the intro- 
duction of a cut-off in the intermolecular potentials 
may well reduce the effective range of the forces, 
producing a different result for the asymptotic be- 
haviour of the unbending line. 

• The intersection of the unbending coexistence line 
with the bulk liquid-vapour coexistence line is tan- 
gential in the BD phase diagram. This feature can- 
not be correct because the difference in coexisting 
adsorptions remains finite at A/^ = 0. 

• The loci of the unbending critical point in the BD 
phase diagram show a nontrivial behaviour as a 
function of the stripe width L\. That behaviour is 
analogous to that of Fig. 7 and is explained above 
in terms of unbending in the presence or absence 
of energetic barrier in the effective potential W\. 



However, the BD phase diagram shows that, addi- 
tionally, for even smaller values of the stripe width, 
the loci of unbending critical points curve towards 
the liquid-vapour coexistence line. To explain that, 
we note that the minimal model prediction for the 
critical temperature, Eq. (|25|), does not depend on 
the thickness of the adsorbed layer on the surround- 
ing substrate while the critical chemical poten- 
tial, Eq. (p6|), does (in fact, in very simple manner). 
The curving of the loci of unbending critical points 
towards the liquid-vapour coexistence line for very 
narrow stripes, found by BD, is therefore a direct 
consequence of the thinning of the adsorbed layer 
on the surrounding substrate. This illustrates the 
merit of first keeping the value of t\ constant be- 
fore allowing for its own dependence on tempera- 
ture and chemical potential. 



IV. PERIODIC STRIPES 

A. General Considerations 

Consider next a periodic array of stripes made of two 
different materials (Fig. 1(b)) whose adsorption proper- 
ties can be described, when pure, by the effective poten- 
tials Wi (£) and W2 (£) , respectively. The widths of the 
stripes are denoted L\ and L2, and the period L = 
This system can be considered as the simplest prototype 
for studying heterogenous wetting since for a single stripe 
the wetting properties are completely determined by the 
embedding material. Thus, the periodic system can be 
used to study the dependence of the wetting temperature 
on different factors such as composition, degree of domain 
separation, etc. This has already been done to a certain 

14| , but this ap- 
ring phenomena 



degree within density functional theory 
proach seems more suitable to study pad 
close to the heterogeneities of the surface. As regards 
the global wetting properties, we believe that a simpler 
approach will suffice. This system has also been studied 
experimentally [Q , mainly in the context of contact an- 
gles and Cassie's empirical law [^9). This law states that 
the contact angle 6 of a macroscopic drop placed on a 
planar heterogeneous substrate satisfies 



cos 8 — 7^ cos 9i 



(28) 



where Ot is the contact angle of a droplet on a (pure) 
material of type i, and 7^ is the fractional area of the 
substrate made of that material. The implication of this 
law for the wetting temperature of a heterogeneous sys- 
tem is straightforward: 



T w = max{TW,T4 2 ),...}, 



(29) 



(i) 

where X"w is the wetting temperature of the substrate i 
(when pure) . Here we check the validity of this law with 
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the help of the minimal model put forward in section II. 
First, however, we recall briefly some general considera- 
tions concerning this system: 

• For any given period L, the limiting cases L\ — > 
and L 2 — * must produce the pure systems 1 and 
2, respectively. This fact follows trivially from the 
Hamiltonian (^J). Observe, however, that it contra- 
dicts clearly Cassie's law since Eq. (^9|) states that 
infinitesimal amounts of a substance can change the 
wetting properties of the embedding substrate. 

• If the thickness of the adsorbed layers on both infi- 
nite substrates is the same {l\ =£2)1 the functional 
(^|) is fully minimised, term by term, by a flat so- 
lution £{x) = l\=^, as in the single stripe system. 
For that reason, the loci of the points for which 
P[ = £ 2 are not to be crossed by any surface phase 
coexistence line. These are homogeneous states on 
a heterogeneous substrate. Note that since both 
£\ and £\ vary with the temperature and chemical 
potential, the loci of the points for which £\ = £ 2 
may not be simple in the T—fi plane. Indeed, the 
line may not exist at all. 

• If an unbending first-order transition takes place 
between different surface phases, the coexistence 
line in the T — fi plane will verify the Clausius- 
Clapeyron-like equations (||) and (Q). 

B. Minimal Model 

As described in section II, the free-energy of this sys- 
tem is given, within the minimal model prescription, by 
the two- variable function H{£\ 1 £ 2 ) 1 Eq. (||). This free- 
energy depends on the average thicknesses £\ and £2, 
whose equilibrium value, for a given temperature and 
chemical potential, is determined by minimisation of the 
free-energy. This yields the following coupled equations: 

2k(£ 1 -£ 2 )+L 1 W{(£ 1 ) = 

(30) 

2 K (£ 2 ~£ 1 ) + L 2 W 2 , (£ 2 ) = 0. 

Adding both equations we obtain the sum rule 

L\ W[(£i) + L 2 W' 2 {£ 2 ) = 0, (31) 

which is the counterpart of another, obtained with the 
full interfacial model (||) 

/ dx W[{£) + I dx W 2 {£) = 0. (32) 

Eqns. (|3(]) can be solved by means of a graphical con- 
struction. The solution corresponds to the intersection 
of the curves 



£2 = £i + ^Wi(£ 1 ) 

(33) 

£i=h + ^W 2 , (£ 2 ), 

in the £r~£ 2 plane, as seen in Fig. 9. As before, if there are 
multiple solutions, the stable one can be discriminated 
upon comparison of their free-energies, Eq. (^|) . 

Among these solutions, it is straightforward to identify 
the spinodal points since they occur when the two curves 
of the graphical construction are tangential to each other 
at a given point. This leads to the following equation 

(i+£ww)=i. (34) 

which can be obtained, alternatively, from the vanish- 
ing of the determinant of the Hessian of the free-energy 
H(£j,£ 2 ). In addition, possible critical points can occur 
if (Hj) is satisfied at a point of inflexion on any of the 
curves: 

W{"(h) = 

or 

W' 2 "{£ 2 ) = 0. 

(35) 

The existence of two alternative conditions for the pres- 
ence of a critical point suggests that, in principle, there 
could be two different critical points in the T—fi phase 
diagram. This is quite different to that obtained for the 
single stripe system. In particular, as the limiting cases 
L 2 = and L\ — are each associated with one of the 
equations (|35|), it follows that there must be two un- 
bending coexistence lines which, in principle, could occur 
simultaneousy in the T—fi phase diagram. All these fea- 
tures, which emerge naturally from only elementary con- 
siderations of the minimal model, would be much more 
difficult to extract from the full interfacial description 




C. Cassie's Law 

Before applying the model to some specific examples, 
we make connection between Cassie's law and the mini- 
mal model. To do that, we define the fractional area of 
material 1 as ^ = L\jL for a given period L. This allows 
us to write the free-energy of the system per unit length 
in the following way 

h{h,h) = j(£i~£ 2 ) 2 + iW 1 (£ 1 ) + (i- 7 )w 2 (4 2 ), 

(36) 

which, in turn, enables us to change from a description in 
terms of L\ and L 2 into another in terms of L and 7. In 
this latter case, L can be understood as a measure of the 
overal size of the stripes, while 7 will be the fractional 
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area of material 1 . This is obviously true as long as 7 is 
not too close to the limiting cases 7 = and 7 = 1. These 
limiting cases give rise (trivially) to the pure substrates 
2 and 1, respectively (for any value of the period L). 

Consider now, that the period L is very large. This 
represents systems for which the heterogeneities extend 
over macroscopic areas. In this case, we can neglect the 
first term of ( |36| ) and, consequently, the variables £\ and 
£2 decouple. Minimisation of the energy is achieved by 
i\ = £\ and £2 = £2 and the free-energy of the system, 
w{T, (i) = vam. h, can be approximated 



u(T,n) « jWx(£l) + (I-7) W 2 (%), 



(37) 



which is but a different way of writing Cassie's law, Eq. 
(p8j) , as generally accepted. The minimal model, there- 
fore, recovers Cassie's phenomenological law in the limit 
of macroscopic heterogeneities. 

The other limiting case, corresponding to small val- 
ues of L, represents systems for which the hetero- 
geneities are of microscopic size. These systems are quasi- 
homogenenous since the two materials are closely inter- 
mingled on the surface. Recall that our model was not 
intended to work under these assumptions. In spite of 
this, the limit L — > produces a sensible result. This 
occurs because the first term of the free-energy (^) be- 
comes exceptionally important and has to be minimized 
fully yielding £\ =£2', i- e. the thickness of the adsorbed 
layer is approximately the same at every point of the 
substrate. The system behaves, in this limit, as a homo- 
geneous system whose adsorption properties will be given 
by the following averaged effective potential 



W(£) = jW 1 (£) + (l-f)W 2 (£). 



(38) 



This description tallies with the spirit of the construction 
of the effective potentials as an integral of the solid-liquid 
and liquid- liquid interactions over the entire system [p0| . 

Therefore, the minimal model seems to describe the be- 
haviour of the heterogenous system in the two mentioned 
limits (macroscopic and microscopic heterogeneities) in a 
physical way. We are fully aware of the fact that such a 
broad description is achieved by smearing out some sub- 
tle details, and that the model will not work under cer- 
tain special circumstances. Nevertheless, we stress that 
the model is put forward just as a versatile approxima- 
tion to the full effective interfacial description, Eq. (fy). 



D. Examples 

As suggested above, we expect a very rich behaviour in 
the periodic system of parallel stripes. A full exploration 
of the phenomenology of this system is a colossal task, 
well beyond the scope of this paper. Therefore, we only 
report here the calculation of the wetting temperature for 
two heterogeneous systems as an example of the method. 



We consider substrates undergoing first-order wetting 
transitions (when pure), whose effective potentials will 
be modelled by 



w^i) = 

and 

W 2 (£) = 



4Bi 



(T — T^) I — 



1 



A, 



£p £p+ 1 



B 1 

£p+- 



(39) 



4B 2 ' 1 w 7 £p 



A-2 



B 2 



(p+2 ' 

(40) 



Out of bulk liquid-vapour coexistence, the usual term 
A[i5p£ must be added, although we do not consider this 
possibility here. The parameters A±, B\, A2, and B2 
are (positive) Hamaker constants characterising the per- 
tinent potentials. Notice that a dimensionless constant 
a has to be included (multiplying the temperature in 
W^ii)) to account for differences in surface entropy be- 
tween the two pure systems. We use A = (Ai/n) 1 / 4 as 
unit of length. 

From this vast space of parameters we have chosen 
two examples, both for long-range (dispersion) forces, 
p = 2. For each of them, we have calculated the wet- 
ting temperature of the periodic system as a function of 
the fractional area of material 1, 7, for different values 
of the period L. As mentioned above, the modification 
of the wetting temperature with the composition of the 
substrate constitute a distinctive and important charac- 
teristic of these periodic systems since, for a single stripe, 
the wetting temperature is necessarily determined by the 
(infinitely wide) abutting substrate. The results, for dif- 
ferent choices of parameters, are shown in Figs. 10 and 
11. 

Let us concentrate first on Fig. 10. As expected, the 

(2) 

wetting temperature ranges between Tw (for 7 = 0) and 

t4 (for 7 = 1), independently of the value of the period 
L. For small periods, the wetting temperature interpo- 
lates almost linearly between Tffl and (as a func- 
tion of 7). This linearity is merely fortuitous. In general, 
T w {l) will be a non- linear function joining the limiting 
cases Tw 2 -* and . Within the minimal model, the wet- 
ting temperature of the heterogeneous system can be cal- 
culated, in this limit L— >0, as the wetting temperature of 
the averaged potential W(£), Eq. (Bq). More specifically, 
for the model effective potentials ( |39| ) and (|4(j), we can 
write an analytical expresion for the wetting temperature 
in the mentioned limit 



Tw(7) 



7 r^+a(i- 7 )r, 

7+a(l — 7) 



(2) 



(41) 



7(l-7)(AiB 2 -A 2 Bi 



(7+a(l-7))BiB 2 (7Bi + (l-7)B 2 )' 
which is clearly non-linear. 
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As the period is increased, Cassie's law, Eq. (|2^), is 
recovered and T w (7) « Tw^ > T w 2 ' , for almost every value 
of the fractional area 7. Deviations occur only for small 
7. These deviations become progressively less important 
for large periods L, as expected. 

The same quantity, calculated for a different system, 
is plotted in Fig. 11. In this case, the differences between 
substrates 1 and 2 are stonger: the thicknesses of the 
adsorbed layers on the pure substrates become very dis- 
similar for temperatures T < Tw and so do their interfa- 
cial tensions (free-energies). In particular, the adsorbed 
layer on substrate 1 is (in general) thicker than the cor- 
responding layer of substrate 2 (see Fig. 12) and its free 
energy is lower. This introduces new phenomenology in 
the phase behaviour. For small periods, L — 0.1 A, the 
wetting temperature of the periodic system resembles an 
inversed Cassie's law since r w (7) « T w 2) <T W 1) , for most 
values of the fractional area 7. The fact that the ad- 
sorption properties of the substrates differ considerably 
appears to have a decisive influence on the adsorption 
properties of the system in the quasi-homogeneous limit, 
L^O. As we increase the value of L, we begin to recover 
the behaviour dictated by Cassie's law. On route to this, 
however, we uncover unexpected richness. For the pa- 
rameters used here, this first occurs for L>0.2A, and is 
clearly visible for L— 1.0A. The broken line in Fig. 11 
corresponds to a (first-order) unbending transition where 
adsorbed layers of different thicknesses coexist (similar to 
those occurring for a single heterogeneous stripe). This 
coexistence finishes at an unbending critical point (rep- 
resented by a black circle) as it did in the single stripe 
system. Of course, in the periodic system, the unbending 
phase transition is associated with an adsorption jump 
on both domains. Fig. 12 shows the thicknesses of the 
adsorbed layers £\ and £2 on substrates 1 and 2, respec- 
tively, as a function of the temperature for a particular 
system (L = 10A, 7 = 0.1) corresponding to those wet- 
ting temperatures plotted in Fig. 11. The thickness of 
the layers adsorbed on the pure substrates l\ and £\ are 
also shown for the sake of comparison. Observe that for 
(T - T^X/At > -0.02 the liquid interface has unbent 
and is essentially flat (£\ w £2). Bearing this in mind, 
we refer to the phase diagrams of Fig. 11 and notice (for 
instance, for L = 1.0A or 10A) that the presence of the un- 
bending transition modifies the qualitative shape of curve 
representing the wetting temperature T w ("/). Specifically, 
at the point where the unbending transition merges with 
the curve T w (j), the latter has a kink. At the right of 
that point, T w (7) grows with a slope similar to the slope 

of the unbending coexistence line, and tends to very 
rapidly (as a function of 7). We can interpret this as an 
unbending-mediated wetting transition. In other words, 
at T w ("f) the interface unbends to the wet configuration 
£\ = £2 — 00. On the other hand, to the left of that 
merging point, the interface finds a stable (unbent) flat- 
ter configuration (£i~£2) resembling those of the quasi- 
homogenenous limit (L«0). Observe that T w (j) hardly 



changes by increasing L in that region of the phase dia- 
gram from the curve obtained for L = 0.1A. This region, 
and the unbending line, are squeezed towards small val- 
ues of the fractional area 7 for larger values of the period 
L. Notice that the unbending line can neither dissappear 
(because that limit contains the unbending transition for 
the single stripe line) nor merge with the line 7 = (repre- 
senting pure substrate 2, which does not undergo any un- 
bending transition). However, it becomes progressively 
and proportionaly less important as L grows and remains 
confined to the region of small fractional area, 7 « 0. In 
the limit of macroscopic heterogeneities, L— >oo, Cassie's 
law is recovered as discussed earlier. 



V. CONCLUSIONS 

In this paper, we have forwarded a simple model of 
surface phase behaviour for fluid adsorption on hetero- 
genenous micropatterned substrates. This has been ap- 
plied to two proto-typical systems: a) a single stripe of a 
material exhibiting either a first- or second-order wetting 
transition embedded in an infinite partially wet substrate 
and b) a periodic array of stripes each of which exhibit 
first-order wetting at distinct temperatures. The advan- 
tage of our " minimal" description is that it allows us to 
explore the vast phase space of parameters that these 
systems present. This would be an extremely arduous 
task within even an effective interfacial Hamiltonian de- 
scription, let alone a microscopic density functional ap- 
proach. Despite the simplicity of our model, we rederive 
results obtained previously using more microscopic meth- 
ods for specific choices of interfacial binding potentials. 
This gives us confidence concerning the validity of our 
method. Applying the minimal model to other choices 
of substrate heterogeneity, we predict new examples of 
surface phase behaviour. In particular, for the single 
stripe system with first-order wetting, we establish the 
conditions under which the unbinding coexistence line 
approaches the pre-wetting line (of the pure substrate) 
from above or below, in the limit of large stripe widths. 
For the periodic system, we have concentrated on estab- 
lishing the value of the wetting transition temperature 
as a function of the fractional area of the two materials 
as well as the total period of the heterogeneity. Agree- 
ment with Cassie's empirical law is obtained in the limit 
of macroscopically large domain sizes. Thus, the wetting 
temperature of the heterogeneous substrate is essentially 
determined by the maximum wetting temperature of each 
of the pure components. For mesoscopic systems, on the 
other hand, notable deviations from Cassie's law are pos- 
sible. This includes "inversions" of Cassie's law where 
we observe that the wetting temperature of the hetero- 
geneous substrate is close to the lower of the two wetting 
temperatures of the pure components. The value of the 
wetting temperature may also be sensitive to the presence 
of unbending phase transitions induced by the heteroge- 
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neous substrate. It is clear to us that even within the 
context of the present minimal description, the possible 
surface phase behaviour on the periodic micropatterned 
substrate is extraordinarily rich. For example, out of 
bulk two phase coexistence, additional phenomenology 
may arise due to the competition between lengthscales 
associated with the pre-wetting line on each pure com- 
ponent. A full exploration of the possible surface phase 
behaviour requires much further work along these lines. 

We are very grateful to Prof. M.E. Cates for support 
during the completion of this project. C.R. acknowl- 
edges economical support from the E.C. under contract 
ERBFMBICT983229. 
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APPENDIX A: 

Consider an infinite stripe of width L\ made of mate- 
rial 1 that is embedded in a substrate made of material 2 
(Fig. 1(a)). If these materials are such that their adsorp- 
tion properties can be described (when pure) by effective 
potentials Wi and W 2 , the adsorption properties of the 
combined system can be obtained by the functional min- 
imisation of the free-energy (0). If substrate 1 undegoes 
a first-order wetting transition at T = when pure 
(and, therefore, its phase diagram presents a pre-wetting 
line), the combined system must undergo an unbending 
transition whose coexistence line must merge with the 
pre-wetting line of the pure substrate in the limit L \ — > 00 
(as explained in the main text of the paper). 

In this appendix, we want to determine under which 
circumstances the interfacial model (||) predicts that, in 
the mentioned limit L\ — ► 00, the unbending line ap- 
proaches the pre-wetting line from above or from below 
(i.e. from temperatures higher or lower than the wet- 
ting temperature , respectively). We denote these 
two possibilities as scenarios A (above) and B (below). 
As we are mainly interested in the merging of the two 
lines close to the bulk liquid- vapour coexistence line (see 
HIE), we only need to concentrate on A/i = 0. 

To start, we notice that our goal can be achieved, with- 
out lack of generality, by studying the system at the wet- 
ting point T=T"4 . Scenario A will be correct if the sta- 
ble configuration, at that point, corresponds to a bound 
state in the limit L\ — > 00. On the other hand, if the 
stable configuration is unbound in that limit, the sce- 
nario B will be followed. In other wor ds, for T = ri 1} 
and A/i = 0, there must be two different configurations 
which minimise the Hamiltonian (^) in the limit L\— >oo 
(see below). This occurs due to the presence of a local 
maxima (the activation barrier) in the potential W\(£) 



at the wetting point. We denote these configurations by 
£~(x) (bound) and^ + (cc) (near-unbound). li£~(x) is the 
stable configuration (it has lower energy), scenario A will 
be correct. Otherwise, it will be scenario B. 

Minimization of (j|) yields the following Euler- 
Lagrange equation: 



W{(£) if \x\ < 
W' 2 {£) if |z| >Li/2, 



(Al) 



where we have located the origin x — in the center of 
the stripe. This equation can be integrated once to yield: 

y ( AWi(£)-V if \x\ < Li/2 

-£(xf = { (A2) 

1 (aw 2 (£) if |x] > £1/2, 

where AWi(£) = W^f) - Wi(£%), AW 2 (£) = W 2 (£) - 
^2(^2)1 an( i is a (non- negative) constant of integration 
to be determined. Note that the boundary condition at 
x = oo has been imposed. Further conditions read 



AWi(£ Q ) = V 
AWi(£ L ) = AW 2 (£l)+V, 



(A3) 



where £0 =£(0) and £l = £(Li/2). These two equations 
can be solved by means of a geometrical construction for 
an arbitrary value of V (see Fig. 13). The relevant value 
for our case (the limit L\ — > 00) is V — > (see below). 
Notice in the figure that the different pairs of values £q , 
£^ and £q , £\ will give rise to the two solutions £~(x) 
and £ + (x). 

Here we can make our first prediction with the help 
of the geometrical construction: if £\ > £ 2 , the correct 
scenario must be A for any physical shape of the effective 
potentials W\ and W 2 . This happens because, under 
these circumstances, every term of the Hamiltonian turns 
out to be lower for the bound configuration £~(x) than 
for the near- unbound one (if the latter exists at all). This 
(purely geometrical) feature is captured correctly by the 
minimal model (^|). However, if £\ < £ 2 at the wetting 
point, the two configurations exist (in this case £ + (x) 
always exists for sufficientely large values of L\) and their 
free-energies need to be calculated. 

To proceed, we use equation ( A2 ) to write the quadra- 
ture 



Li 



dx 



y/AW!(e)-P 



(A4) 



which determines V implicitly as a function of the stripe 
width Li, T and /i. A simple calculation shows that V is 
related to the free-energy of the system in the following 
way: 



V 



9Z7 



(A5) 
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Thus, the comparison of the energies of configurations 
£~{x) and £ + {x) can be obtained from the asymptotic 
behaviour of V for large values of L\. In turn, t his follows 
from studying the asymptotic behaviour of (A4) when 
V — >0, since this is the only value for which the integral 
diverges and recovers, in that way, the limit L\ — > oo. 
For appropriate effective potentials W\ and W% (with the 
only requirement that W\ has an activation barrier), we 
obtain 



2p 



(A6) 



where (— ) and (+) denote the bound and near-unbound 
configuration, respectively, and p is the range of the 
forces, i.e. W\{1) <~ £~ p when I — ► oo. Therefore, the 
energy of each configuration can be written asymptoti- 
cally as 



3-P 



(A7) 



where u>q , u>q, cjF and wj~ are constants that depend 
on the details of the potential. Note that the irrelevant 
common term [Wa^) — Wi. (■££)] L\ has been subtracted. 

Consequently, the range of the forces p appears to be 
determinant in establishing which of the two scenarios 
holds. For p < 3, the energy of the bound configuration 
tends to a constant for large stripe widths, lj~ — > loq , 
while the energy of the near-unbound configuration di- 
verges as a function of the stripe width L%. For that 
reason, the bound configuration £~(x) is the stable one 
in this case, and scenario A is correct. This includes the 
experimentally relevant case of dispersion forces, p = 2. 
In contrast, for p > 3, both w~ and u> + tend to constant 
values in the limit L\ — > oo. The prevalence of scenario 
A or B will be determined, in this latter case, by the 
comparison of the free-energies: 



2V2E 



d£ \J AW\ {£) + / d£ y/AW 2 (e), (A8) 



and 



2V2E 

The necessary conditions for scenarios A or B to occur 
are gathered in the following table: 



Conditions 


Scenario 




A 




p<3 
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p>3 
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UJ„>UJ+ 


B 



where P[ and £\ are the wetting layer thicknesses on 
the pure substrates measured at T = T^P . 
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Material 1 Material 2 





Li L 2 



FIG. 1. Schematic illustration of micropatterned sub- 
strates considered in this paper. The figure shows a) a single 
stripe of width L\ embedded in substrate 2 and b) a periodic 
array of stripes of widths L\ and L2, made of materials 1 
and 2, respectively. Both systems are translationally invari- 
ant along the y direction. 



Critical Point 




-\ 


A, Inflexion 







Spinoidal 
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FIG. 2. Sketch of the graphical construction used to solve 
Eq. (JTT]) . The average thickness of the adsorbed layer on the 
stripe, £1, is obtained by the intersection of the curve W{(£) 
with the straight line 2k(P{ — l)/Li (see (a)). The geometrical 
conditions for a critical point (b), a spinodal (c), and the 
coexistence between two different values of the adsorption l\ 
(d) are also shown. 
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FIG. 3. T—fj, phase diagram for the adsorption on a stripe 
of different widths L\. The stripe is made of a material 
that, when pure, undergoes a second-order wetting transi- 
tion at T — Tw 1 ' (large full circle) , whose adsorption proper- 
ties are described by the effective potential (|14[). The stripe 
is embedded in a partially wett substrate with finite film 
thickness £J = 4A. This system shows an unbending tran- 
sition (black continuous lines) for a range of stripe widths. 
In the figure, Li/X = 12,16,20,24,28,32,36,40,44 and 48. 
For L\ » 53 A, the unbending transition is pre-empted by 
the bulk liquid- vapour phase transition (large empty circle). 
The continuous grey line represents the loci of the critical 
points (as given by Eqs. ( |l5[ - |l^ ). This line extends into the 
metastable region of the phase diagram (dashed grey line) and 
tends asymptotically to the wetting temperature T,t when 
L\ — » 00. Note that the unbending transition does not cross 
the iso-adsorption line l\ — £2, (dotted line). 
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A^Sp B 4 /A 5 



FIG. 4. pre- wetting line (black) as modelled by the effective 
potential (BOj) withp = 2. Iso-adsorption lines (grey) are also 
shown for different values of the adsorption £ ranging between 
1.9 and 20 (in units of B/A). 
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FIG. 5. T—/J, phase diagram for the adsorption on a stripe 
of width L\ made of a material exhibiting a first-order wet- 
ting transition (when pure). The stripe is embedded in a 
partially wet substrate whose wetting layer thickness t\ is 
considered independent of temperature and chemical poten- 
tial. The pre-wetting line of the (pure) material 1 is shown 
for comparison (See text for details). 
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FIG. 6. T— p phase diagram showing unbending phase tran- 
sitions occuring on a substrate with a single stripe showing 
a first-order wetting transition (when pure), embedded in a 
substrate whose wetting layer thickness = 10A. The un- 
bending coexistence lines are represented by thin black lines 
which finish at critical points, the loci of which is shown as 
grey. Different stripe widths are considered: L\ = 1,2,5, 10 
and 20 (xl0 3 A). The pre-wetting line of the (pure) mate- 
rial 1, which corresponds to the unbending line in the limit 
L\ — > oo, is shown for comparison (thick line). The broken 
line represents the line for which P[ =1%, Observe that this 
line is not crossed by any unbending line (See text for details). 
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FIG. 7. As Fig. 6, with £J = 4A and U = 2, 3, 4, 5, 7.5, 10, 20 
and 50 (xl0 2 A). Inset: Merging of the unbending line with 
the pre-wetting line for large stripe widths, L\ — 1, 2, 4, 8 and 
16 (xl0 4 A). 
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FIG. 8. As Fig. 6, with 11 = 10A/3 and different values of 
L\. The grey broken line is the loci of the unbending critical 
points that have been pre-empted by the bulk liquid-vapour 
phase transition. Inset: Merging of the unbending line with 
the pre-wetting line for large stripe widths, L\ = 4, 8, 16, 32 
and 64 (xl0 3 A). 



14 









~/\* I 1 = I 2 +(L 2 /2k)W 2 



FIG. 9. Sketch of the graphical construction used to solve 
equations In the figure, both potentials W\ and Wi 

represent first order wetting transitions. The solutions corre- 
spond to the intersection of the curves £2 — £i + j^W{(£i) and 
£\ =^2+2^2 (^2) (black circles). If there are more than one, 
as in the figure, the stable one will have the lower free energy, 
Eq. (|^). The thickness of the adsorbed layers on the pure 
substrates £\ and £ 2 correspond to a crossing of the curves 
with the line £±=£2 (white circles). 
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FIG. 10. Wetting temperature (measured with respect 
to ri, 2 ') of a periodic heterogeneous substrate calculated 
with the minimal model, Eq. (||) and the effective poten- 
tials ( p9| ) and (^). The wetting temperature (in units of 
Ai/X) is represented as a function of the fractional area 
of the first material, 7, for different values of the pe- 



riod L/X = 1,10, 10 2 , 10- 



,10 4 ,10 5 , 



10 6 and 10 7 (from be- 



low). We use the following parameters: T 5 
B 1 = 2B 2 =A 1 X, Ax=A 2 and a = 2. 
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FIG. 11. As Fig. 
B 2 = ABi = 2AiX, A 2 = 0.05 Ai and a = 20. The broken 
fine represents an unbending coexistence line appearing for 
L~0.2A (see text for details). 
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FIG. 12. Adsorption on the periodic system whose wetting 
temperature is plotted in Fig. 11 for L= 10X and 7 = 0.1. The 
curves £\ and £2 represent the adsorption on the substrates 
1 and 2 (respectively). A (first-order) unbending transition 
takes place for (T - T^)X/Ai ~ -0.02. For the sake of com- 
parison, the adsorption on the pure substrates 1 and 2 is also 
shown (£\ and 11). 
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FIG. 13. Geometrical construction to solve Eq. (A3) for 
arbitrary effective potentials Wi and Wi at T = Tjp~- Wi 
shows an activation barrier, necessary to model a first-order 
wetting transition. There are two possible solutions for both 
lo (empty circles) and £l (black circles) . These correspond to 
bound (-) and near-unbound (+) interfacial configurations. 
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